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Abstract. In the last fifteen years great progress has been made in the under- 
standing of nonlinear resonance dynamics of water waves. Notions of scale- and 
angle-resonances have been introduced, new type of energy cascade due to non- 
linear resonances in the gravity water waves has been discovered, conception 
of a resonance cluster has been much and successfully employed, a novel model 
of laminated wave turbulence has been developed, etc. etc. Two milestones 
in this area of research have to be mentioned; a) development of the g-class 
method which is effective for computing integer points on resonance manifolds, 
and b) construction of marked planar graphs, instead of classical resonance 
curves, representing simultaneously all resonance clusters in a finite spectral 
domain, together with their dynamical systems. Among them, new integrable 
dynamical systems have been found that can be used for explaining numerical 
and laboratory results. The aim of this paper is to give a brief overview of our 
current knowledge about nonlinear resonances among water waves, and finally 
to formulate the three most important open problems. 

1. Exposition. In this paper we wiU try to present a major part of known analyt- 
ical, numerical and laboratory results on nonlinear resonances among water waves, 
in as strict mathematical language as possible. This is not a simple task due to the 
three-fold problem: 1) there is no strict definition of a wave; 2) there is no general 
agreement about the types of waves which should be called water waves; 3) the 
notions of resonance in physics and mathematics are different. Let us go through 
all these points one by one, regarding for concreteness 2D-wavevectors. 

First, the simplest possible understanding of a (propagating) wave as a Fourier 
harmonics 

exp z(k • X — cij t) (1) 

is obviously too simplified and does not include normal modes which are due to 
boundary conditions. Here x = (xi, X2 ) and time t are space and time variables 
correspondingly, to = w(k) is dispersion function and k is wavevector. For instance, 
the normal mode of oceanic planetary waves (that are due to the Earth rotation) 
with zero boundary conditions in a finite box [0,Lj:J x [OjLa;^] reads [43] 

/ mxi \ , ( nx2 \ . f P \ 
\A]^\svtiyK— — J sm|^7r-^- — j svay—xi+ujtj. (2) 
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Boundary conditions on a sphere or a circle lead to even more complicated forms 
of the normal mode including special functions [43]. "We seem to be left at present 
with the looser idea that whenever oscillations in space are coupled with oscillation 
in time through a dispersion relation, we expect the typical effects of dispersive 
waves" [61]. So, from now on wc assume that dispersion function a;(k) defines 
the type of a wave, for instance, lo ^ |k|^/^ corresponds to capillary waves in 
a rectangular box with periodic boundary conditions, lo ^ |k|~^ corresponds to 
oceanic planetary waves in a rectangular box with zero boundary conditions, etc. 
Sign "~" instead of "=" in definitions of dispersion function means that some 
constant is omitted. 

Second, capillary, gravity-capillary, gravity surface, oceanic planetary and freak 
waves, tsunami, internal waves in rotating fluid, etc., etc. can be observed in water, 
and they all can be - and sometimes are - called water waves. This is, of course, 
too extensive a definition that does not allow any reasonable overview, neither in 
the form of a paper nor of a book. An example of another extreme can be given by 
[37] where only waves on the surface of deep water are called "water waves" . In this 
paper we choose to follow the classical conception of water waves [61], as defined 
by dispersion function 

= (.gfctanhfcd)(l + — fc2), with A: = |k| (3) 

and g, d, (7, p being the gravitational constant of acceleration, the average water 
depth, the coefficient of surface tension and density correspondingly. In this paper 
we will concentrate on, though not restrict ourselves completely by, the three lim- 
iting cases of (3) for water waves, with p = 1 and fed — > oo, which are important 
for numerous applications: 

1. capillary waves, Xch 2 cm : oj^ = ak^, (4) 

2. gravity waves, Xch S> 2 cm : lo'^ = gk, (5) 

3. gravity-capillary waves, Xch ^ 2 cm : lo'^ — gk + crfc'^, (6) 

where Xch denotes characteristic wave length for each approximation. 

Third, the notion of resonance appears in mathematics in the theory of Poincare 
normal forms, as a condition of linearization of a nonlinear ODE or a system of 
nonlinear ODEs [42], in the form 

PIUJI H ^Psi^s +Ps+lUJs+l = 

where loj is notation for w(kj) and pi, p2, ■ ■ ■ , Ps+i G ^ are integers. The physical 
notion of resonance has been first introduced by Galileo Galilei [16] who studied 
oscillations of a pendulum under the action of a small external force. In modern 
language, the corresponding equation reads xn + r^x = ecxp{ift) with < e <C 
1, and the oscillation amplitude grows linearly with time if eigenfrequcncy of the 
system lo coincides with the frequency of the driving force /. The analogy with the 
three-wave resonance is now obvious: two Fourier harmonics of the form (1) with 
wavevectors ki, k2, frequencies loi, lo2 and amplitudes Ai, A2 will give resonance 
with harmonic 

A3 exp [i(ki ± k2) • X - {uji ± UJ2) t], 
which yields 3- wave resonance conditions in the form 



L0i+L02=L03, ki+k2=k3. 



(7) 
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Similarly, 4- wave resonance conditions read 

UJi + U2 = ^3 + ^^4, ki+k2 = k3 + k4. (8) 

In the case of periodic or zero boundary conditions, solutions of (7), (8) have to be 
found in integers, i.e. for kj = {mj,nj) with mj, nj S Z. Solutions of (7) and (8) 
are called triads and quartets correspondingly, they are defined by the geometry of 
the wave system and do not depend on time, that is, they describe kinematics of 
a wave system. Dynamics of a resonant triad is covered by 

iA, = V^l*A;A3, zi2 = ViV^J^3, iAs = V,\A,A2, (9) 

where Aj are amplitudes of three waves satisfying (7) and A = -^A. Dynamical 
system for a resonant quartet reads 

i ii - T^^A^AsAi + [Qi - LOi)Ai , i = T^iAlAsA^ + {^2 - ^^2)^2 , 

i A3 = TifAlAiA2 + {^3 - ^3)^3 , iAi = Tii*AlAiA2 + (n^ - ^^4)^4 . (10) 

where fij are so-called " Stocks-corrected" frequencies [53] . They are often omitted 
in the literature while suitable renormalization of iv-s yields more usual form of 
(10), without terms {flj — ujj)Aj included (e.g. [38], for one-dimensional quartets). 
The general form of (9), (10) can be deduced in the frame of the Hamiltonian 
formalism, [66], and does not depend on the type of waves. The only difference 
between various wave types is hidden in the form of coupling coefficients V12 and 
T3I; their expressions for water waves are given in the Table 1. Resonant triads 
and quartets may form resonance clusters, for instance when one resonant mode 
(to, n) is part of a few solutions of (7) or (8). Obviously, triads and quartets are the 
minimal resonance clusters in three- and four-wave systems correspondingly; they 
are called primary clusters. 

Remark 1. To use the Hamiltonian formalism for describing these wave systems, 
some small parameter < e <JC 1 has to be introduced, which allows the application 
of variational or multi-scale methods. These yield, in the the simplest case, two 
time-scales - fast time, t, for the linear Fourier harmonic (1), and slow time, T = et, 
for the amplitudes of the nonlinear resonant mode (9). In water wave systems, wave 
steepness is usually taken as a small parameter, e ^ 0.1. 

2. Kinematics. 

2.1. Exact resonances. Computing integer solutions of resonance conditions (7), 
(8) for dispersion functions (4)-(5) is equivalent to finding integer points on a res- 
onant manifold of a high degree in many variables. Loosely speaking, the problem 
is equivalent to Hilbert's Tenth Problem [21], which is proven to be unsolvable in 
the general case [40], and even for the case of two variables only partial results are 
known. We are interested in solutions of equations (7), (8) with 6 and 8 variables 
correspondingly. One could hope that numerical search for solutions might be a 
reasonable possibility but it is not. Full search for multivariate problems in integers 
consumes exponentially more time with the size of the domain to be explored and 
for each additional variable. For (8) in the domain |to|, |n| ^ 1000 this would im- 
ply some 10^^ tries. Still worse, the equations include radicals and straightforward 
transformation to a purely integer form would lead to operations with huge integer 
numbers - for the said domain, of the order of lO^^*'. All these reasons make a quest 
for effective algorithms unavoidable. 
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A special method for solving systems of the form (7), (8) has been presented in 
[26] for gravity waves with dispersion function (5). Details of its implementation 
and numerical results for various irrational dispersion functions can be found in 
[29, 30]. The method is called g-class method or g-class decomposition. Its idea 
comes from the linear independence of a finite number of some algebraic numbers. 
Below we give the general definition of the (/-class and show how to use it for finding 
solutions of (7), (8). 

Definition 2.1. For a given c G N, c 7^ 0, 1, — 1 consider the set of algebraic numbers 
Rc = fc^/^, fc G N. Any such number kc has a unique representation 

fc, = 7gl/^7eZ with q = pfpl\..pl\ 

and pi, ...pn being all different primes and the powers ei, ...e„ G N are all smaller 
than c. The set of numbers from Rc having the same q is called 5-class Clq. The 
number q is called class index. 

Now we can decompose the computational domain into g-classes and search for 
solutions in each class separately. The computational advantage of using g-class 
decomposition is huge. For instance, for dispersion function (5) and resonance 
conditions (8) straightforward computation takes 3 days with Pentium-4 in the 
spectral domain |m|, |n| < 128; the g-class decomposition produces all solutions for 
|m|, |n| < 1000 in about 15 minutes with Pcntium-3. 

2.1.1. Capillary waves, uP' ~ ak^ , three-wave resonances. In this case, we take c = 2 
in the definition above and represent the norm kj as 

^3 ^ V^ii 7j' S N, qj is square- free (11) 
The presentation (11) is unique for each kj and it follows that 

uji+uj2^i^z ^ 7i 91 +72 92^92 = 7393 ^/g3■ (12) 

The necessary condition of the existence of integer solutions of (12) is qi ~ q2 ~ qs, 
say qj = q, \/j ~ 1,2,3, (see [25]), and this allows to reduce (12) to a particular 
case of Fermat's Last Theorem, 7^ -I- 7I = 7I, i.e. there are no resonances among 
capillary waves with dispersion function (4). 

The case of capillary waves is unique. Usually, a three-wave system possesses a 
number of resonances divided into small resonance clusters [24] . Mostly, resonance 
clusters are isolated triads or groups of two connected triads (70 — 100%% of all 
clusters, depending on the size of the spectral domain). For instance, the same 
g-class construction used for ocean planetary waves with dispersion function uj ~ 
k~^P and resonance conditions (7), turns the first equation of (7) into 72 73+71 73 = 
71 72, in which case solutions do exist [29]. 

At the end of this section we would like to mention quite recent progress on the 
resonances of capillary waves reported in [8] . Dispersion function has been taken in 
the form 

^ = ^ + ^Vr!2 + 4fc3a, (13) 

where f2 is non-zero constant vorticity. The main results are as follows: 1) wave 
trains in flows with constant non-zero vorticity are possible only for two-dimensional 
flows; 2) only positive constant vorticities can trigger the appearance of three- 
wave resonances; 3) the number of positive constant vorticities which do trigger a 
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resonance is countable; 4) the magnitude of positive constant vorticity triggering a 
resonance can not be too small. Of course, these results are only valid for waves 
with small amplitude [8, 9, 10, 60] (cf. Remark 1 above). 

2.1.2. Gravity- capillary waves, uj'^ = gk+ak^, three-wave resonances. Construction 
of q-classcs is this case differs from (11) and is more intricate [32], due to the fact 
that dispersion function contains two irrationalities and is not homogenous on fc, 
i.e. coefficients g and a do not disappear. For numerical simulations, we used 
g ~ 981 and two values of ct: a ~ 74 and <j ~ 75, for water with temperature 18°C 
and 8°C correspondingly. Numerical simulations have been performed, seeking 
solutions uji + UJ2 = ^3(1 + s), for various £, < |e| < 10~^. Spectral domain 
— 100 < m, n < 100 has been investigated. For both cases, cr = 74 and a = 75, only 
isolated triads have been found (for e = 0): 24 and 16 triads correspondingly; no 
triad appears simultaneously in both lists. The resonance structure is much richer 
for e > 0, resonance clusters of two to six triads appearing already for e = 10^'', 
while for e = 10~^ the overall number of resonances is about 10^, and more than 
90% of all resonant triads coincide for cr = 74 and cr = 75. 

2.1.3. Gravity waves, lo^ = gk, four-wave resonances. In this case, c = 4 in the 
definition of the q-class, i.e. k = {m? + n^)^/^ = 79^^"', where q does not contain 
fourth degrees in its representation as a product of different primes in corresponding 
powers [26]. It can be proven that only two cases are possible: 1) all 4 wavevectors, 
forming a solution of (8), belong to one class Clq^ and 2) two wavevectors belong 
to Clq-^ and two others belong to Clq.^ with qi ^ qi. Again, this is only a necessary, 
not a sufhcient condition for the existence of a solution. 

From physical point of view, the main difference between nonlinear resonances 
in three- and four-wave systems can be formulated as follows. Any three-wave 
resonance generates a new frequency and, correspondingly, new wave length, i.e. 
any three-wave resonance contributes to the energy transfer over the scales. This 
is not true in four-wave systems. Indeed, let us regard one specific solution of 
(8): ki = (-64,-16), ka = (4,16), kg = (-64,16), k4 = (4,-16). This solution 
is not trivial - all wavevectors are different, but k\ = kj, and k^ = ^4, i.e. no new 
wave length appears due to this kind of resonances. Being regarded in k-space, 
these resonances perform the energy transfer not over the scales but over the angle, 
forming two circles. Correspondingly, two types of resonances arc singled out, with 
different dynamics - angle- and scale-resonance [28]. 

The fact of utmost importance is that these two types of resonances should not be 
studied separately, because mixed cascades also exist. For instance, the wavevector 
(119,120) takes part in one scale-resonance and 12 angle-resonance [28]. Complete 
study of all the resonances in the spectral domain |mi|, \ni\ < 1000 can be found 
in [35]. In this domain, there are more than 600.000.000 exact resonances, among 
them only 230.464 angle-resonances, among which 213.760 resonances are formed 
by four coUinear wavevectors. As it is shown in [14], the coupling coefficient (see 
Table 1) in this case is equal to zero, i.e. these resonances do not affect dynamics. 
Some parametric series of solutions for the resonance conditions (8) are known, for 
instance, for any given wavevector (m, n), a 3-parametric series of angle-resonances 
is known: 



ki = (m,n), k2 = (i, -n), kg = (m, -n), k4 = (i,7i), t = 0, ±1, ±2, . . . . (14) 
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Series of scale-resonances called tridents has been first found in [39] for quartets of 
the form: 

ki = (a, 0), k2 = i-b, 0), k3 = (c, d), k4 - (c, -d) , (15) 

where {s'^ + t'^ + stf, {s^ + t'^ - stf, c = 2st{s'^ + i^), d ^ s"^ ~t^ , s,t e N. 
This scries is of special interest for studying large-scale gravity waves: all scale- 
resonances in the spectral domain, say, k < 100 are of this form [35] and therefore 
can be found analytically. Analytical series are very helpful not only for computing 
resonance quartets and cluster structure but also for investigating the asymptotic 
behavior of coupling coefficients . More details can be found in [35]. 

Remark 2. The results reported in this section were possible to obtain because all 
three dispersion functions (4)-(5) are irrational and we could make use of the g-class 
method. Interestingly enough, similar idea can also be worked out into a compu- 
tational method for a transcendental dispersion function, say, = g k tanh k d. 
Indeed, using standard presentation tanh a; = [exp a; — exp (— x)]/[expa;-|- exp (— x)], 
one can rewrite (8) as a combination of different exponents with polynomial coef- 
ficients depending on k and d (see [23], Discussion, p. 52). Afterwards appropriate 
theorems about linear independence of exponents over algebraic numbers can be 
used [2]. This procedure can be used for any transcendental dispersion function 
which has a representation as a rational function of different exponents. 

2.2. Topological structure of the complete resonance set. In 1950-1960th 

the usual way to represent nonlinear resonances was to construct a resonance curve 
that is the locus of pairs of wavevectors interacting resonantly with a given wavevec- 
tor. A locus might be an ellipse, be shaped like an hour-glass, can degenerate into 
a pair of lines, etc. Eventually, the characteristic form of the locus is useful to know 
before planning some laboratory experiments: in [45] a resonance curve has been 
constructed for a special type of resonant quartets of gravity waves, in which two 
wavevectors coincide. A typical resonance curve for gravity-capillary waves can be 
found in [46]. The main drawback of the resonance curve presentation is that it 
allows to visualize - at most - only a part of the resonance set. And of course, 
without the q-class method, there has been no constructive way to establish whether 
any integer points exist on these curves. 

A novel graphical presentation of 2D-wave resonances as a planar graph, [34], 
provides a very clear and transparent way to exhibit the complete resonance set. 
The construction is performed in two steps. First, each 2D-vector is regarded as a 
node of integer lattice in the spectral space, and nodes which construct one solution 
(triad or quartet) are connected. The result is called geometrical structure of the 
resonance set and can be rather nebulous [34]. Second, all different topologically 
equivalent components of the graph, corresponding to the geometrical structure, 
are singled out. The list of all these elements together with a number, showing 
how many times a specific component is met in the solution set, is called topological 
structure of the resonance set. For instance, in a three-wave resonance system, most 
part of the elements are either isolated triads, shown as triangles, or butterflies, 
shown as two triangles with one joint vertex. For dynamical reasons which are 
explained in the Sec. 3.2, it is important to known whether or not a connection 
within a cluster is realized via a mode with maximal frequency (wa-mode). This 
mode is marked by the letter A in each triangle, while two other modes are marked 
by P. Examples of topological elements appearing in the various water wave systems 
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are shown in Fig.l in Sec. 3.2. This construction can be extended to the case of 
four wave resonances but the structure of clusters is more involved [35] . 

3. Dynamics. 

3.1. Primary clusters. A more compact form of (9) can be obtained by a suitable 
change of variables and reads 

Bi = ZB;_B3, B2^ZBIB3, B3 = -ZB,B2. (16) 

It has three conservation laws (e.g. [4]) 

hs = + /i3 = + |B3|2, Htr = lm{BiB2Bl). (17) 

The first two of them, I23 and /13, can be rewritten via energy and enstrophy, 
while the last one Htr is the Hamiltonian of the triad. Notice that though (16) 
has 3 complex variables, i.e. 6 real variables, and only three conservation laws, it 
is easy to show that it is integrable. Indeed, one should rewrite it in the standard 
amplitude-phase representation Bj = Cj exp(i0j): 

Ci = ZC2C3 cos if, C2 = ZC1C3 cos if, C3 = -ZC1C2 cos if, (18) 

(p ^-ZHtr{C^^ +C2^ -C^^). (19) 

where Lp ~ O1+ 62 — O3 \s called dynamical phase. Now the four equations (18), (19) 
have four real variables, and analytical expressions for amplitudes Ci, C2, C3 are 
found in Jacobian elliptic functions dn, sn, en [62]. Using relations between squares 
of elliptic functions, we can represent squared amplitudes as 

CUT) = (^)%n^ [2K{,) +71 

CUT) = (^)'sn2 (2X(/.) n:^,^) (20) 

CUT)^ ^(H^)%n^(2A-(^)(^,^)+73, 

where K{fi) is elliptic integral of the first order, and fi, to, t, 71, 72, 73 are explicit 
rational expressions of the initial values of the amplitudes C^(0), C|(0), C|(0) (see 
[6]). Since the energy of each mode Ej is proportional to its squared amplitude 
Cj and < sn^ < 1, the characteristic energy variation of any resonant mode 
Emode ^ sn^{k t, fi) is bounded: < E^node < Eq, with Eq = £'o(/i3, ^23, Htr) being 
an explicit function of the initial conditions. 

The dynamical phase satisfies an evolution equation (19) and its solution cannot 
be obtained by simply replacing the solution for the amplitudes in the Hamiltonian 
Htr ~ C1C2C3 sintp and solving for cp. The reason is that non-zero ip generically 
evolves between and tt, crossing the value = 7r/2 periodically. This implies 
that sin~ is double-valued and thus it is not possible to obtain 1^9 in a unique way. 
Importance of the dynamical phase is due to the fact that it effects substantially 
the magnitudes of Ci,C2,C3 (see Fig.l from [5]). Solution for ip, first obtained in 
[6], reads 

/ M f'^KU)Y ( , AT -to) 
ip{t) = sign (v3o) ar coot I ^j^^ \ Z~r~ ) ^^^^^^ ^iP-) ~ 

' (21) 

with notation sn en dn(a;, n) = sn{x, fj) cn(a;, fj) dn(a;, /j,) . 

Substantially less is known about the integrability of (10). Its analysis can prob- 
ably be carried out along the same lines as for (9). A promising start has been made 
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in [53] , though analysis was not brought all the way to explicit formulae similar to 
those obtained for (9). Also the importance of the dynamical phase has not been 
put to work yet, while (10) is analyzed in the complex form and the effect of the dy- 
namical phase is hidden. For the general form of a quartet and for arbitrary initial 
conditions, (10) "docs not exhibit strict periodicity" [53], though in the particular 
case of a trident, its behavior is periodic [52] (both statements refer to the results 
of numerical simulations). 

3.2. Clusters of two and more triads. Dynamics of a cluster consisting of a 
few connected triads depends on whether or not the connecting mode is the mode 
with maximal frequency uj^ in one or both triads [33]. It follows from the form of 
conserved quantities /23,^i3 (see (17)) that if only one of the modes with lower fre- 
quencies oji and UJ2 is exited, it can not share its energy with other modes. Therefore 
these two modes are called P-modes, P from passive. On the contrary, wa-mode is 
called A-mode, A from active, because being exited initially, it causes exponential 
growth of P-modes amplitudes, until all the modes will have comparable magnitudes 
of amplitudes (cf. the criterium of instability in [20]). In Fig. 1 examples of simple 
resonance clusters are given, with marked vertexes. Each isolated marked graph 
defines uniquely some dynamical system. For instance, a PP-butterfly consisting of 
triad a and triad 5, connected via the mode Bia = -Bib, is covered by 

Bla = ZaB2aBia + ZhB^ifB^b, B2a = ZaB^^Bsa, 
B2b = ZbBl^B^b, B^a ~ —ZaBiaB2a, B^h = —Zi,BiaB2b ■ (22) 

and conservations laws have the form 

hsa = |-B2a|^ + l-Bsal^ , /23b = |-B2b|'^ -|- |i?3b|^ , 
lab = \Bl\'^ + |-B3a|^ + \B3b\'^ , Hbut — ^^{ZaBiB2aB^g^ + ZbBiB2bB^ij) . 




Figure 1. Triad with marked vertices; PP-butterfly; AP-PP-kite, 
AA-ray and AP-PP-chain of three triads. 

Some scattered results on the integrability of resonance clusters can be found 
in [4, 6, 57, 58]. For instance, PP-butterfly (22) is known to be integrable if 1) 
Za/Zh = 1 or 2 or 1/2; and 2) Za,Zi, are arbitrary and Htut — 0. However, it is 
not proven that this list is exhaustive. Some of the kites, rays and finite chains of 
triads (see Fig.l) are integrable, but complete classification of integrable resonance 
clusters is still to be constructed. 

No quantitative results are known about dynamics of resonance clusters in four- 
wave systems. Simple qualitative scenario suggested in [35] is based on the pure 
kinematic considerations - reservoirs formed by many angle quartets in quasi- 
thermal equilibrium with sparse links between them formed by scale-resonances 
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VI2 

[64] 



87rV2cr 



^WlW2W3 X 



where /Cfe^^fca = (k2 • ks) + fc2fc3 • 



(g/fc + afc3)i/2 



[41] 



-(k2-k3)(t-t-^^) 



5fc 



1/2 



^34 

[63] 



12kik2 k^ki 



-2 
-2 
-2 



[^3^4(ki • ka - fcifc2) +tJiCJ2(k3 ■ k4 - fc3fc4)] 
[cJ2t^4(ki • ks + fcifcs) + CJiLJ3(k2 • k4 + ^2^4)] 

[w2W3(ki • k4 + kiki) + cjicj4(k2 • ka + fc2fc3)] 





a' 














9' 



At,, 1 ,, \2 (ki-k2-feife2)(-k2-k4+*;2fc4) 



+A(UJ1 - t^3)2(ki-k3+fcifc3)(k2-k4+fc2fc4) 



, A(, , , , \2 (ki-k4 + fcifc4)(k2-k3 + fc2fc3) 

+4(a;i-LJ4j 

+ (ki • k2 + fcifc2)(k3 • k4 + fc3fc4) 
+ (-ki • kg + fcifc3)(-k2 • k4 + fc2fc4) 

+ (-ki • k4 + fcifc4)(-k2 • k3 + fc2fc3) 



Table 1. Overall data on three types of water waves: form of 
dispersion function, notation for coupling coefficient and original 
reference, form of coupling coefficient (1st, 2d and 3d columns cor- 
respondingly) . 



- and has yet to be justified. At least, magnitudes of coupling coefficients have to 
be computed; this qualitative scenario will be invalidated, for instance, if most of 
angle-resonances have zero coupling coefficients, or their magnitudes are substan- 
tially smaller then those of scale-resonances, etc. 
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Remark 3. Three- wave system might possess a four- wave resonance cluster but it 
will not be a primary cluster (sec Fig.l, kite or ray). This can easily be seen from 
the form of the corresponding dynamical system. For instance, a PP-PP kite, with 
Bla = Bib and B2a = ^26, has dynamical system 

Bla ~ B2a{ZaBza + Z^B^b) , B2a = B^aiZaB^a + ZbB^b) , 

Bsa = ~ZaBiaB2a , B^b = — ZbBiaB2a , (23) 

which is different from (10) . Clusters of this form are indeed observed among oceanic 
planetary waves [36] and their integrability is proven in [6] . 

Remark 4. All the results above have been obtained for water waves without 
vorticity. As has been shown in [8], non-zero vorticity can generate non- linear 
resonances. It would be interesting to derive dynamical equations for resonance 
clusters in this case. The presence of non-zero vorticity invalidates the existence of a 
velocity potential for the flow and harmonic function theory is not readily available 
for the analysis. On the other hand, both multi-scale methods and variational 
techniques can also be used for the case of non-zero vorticity [9, 11, 17, 59]. 

4. Numerical simulations and laboratory experiments. In the past, numer- 
ical simulations [1, 39, 48, 65] and laboratory experiments [13, 15] have been mostly 
performed for checking the prediction possibilities of classical wave turbulence the- 
ory with continuous wave spectra (CWT), that is, kinetic equations and power 
energy spectra. It has been established that the so-called infinite-box limit as- 
sumed by CWT theory is not achieved, and even in a 10m x 6m laboratory flume 
finite-size effects are very strong [13]. 

Very interesting results are presented in [7], where chains of three connected 
quasi-resonant triads of gravity-capillary waves have been identified in the exper- 
imental data, dynamical system consisting of 7 nonlinear DDEs was written out 
explicitly and solved numerically. Experimental data have been compared with the 
numerical predictions of 7-modes model computed for constant dynamical phase, 
(y9 = TT. Predictions were in qualitative but not in quantitative agreement with the 
experiment (magnitudes of amplitudes were underestimated by the theoretical pre- 
diction). Quite recently results of laboratory experiments with surface waves on 
deep water were reported [22, 50] in which regular, nearly permanent patterns on 
the water surface have been observed. It would be interesting to attribute these 
regular patterns to specific resonance clusters, like it has been done in [7] but taking 
into account the dynamical phase also. 

In this context, we would like to mention results presented in [49]. Pattern 
formation due to nonlinear resonances is studied numerically for a somewhat sim- 
plified model PDE (nonlinear terms do not include any derivatives). Dissipation 
and multifrequency forcing are taken into account; the importance of the modes' 
phases is fully recognized; pattern-forming modes are computed via fast Fourier 
transform. The main drawback of this Ansatz is that the model equation, though 
it has a Hamiltonian limit, can not be derived from Navier-Stokes equation and its 
dispersion function is different from the real-world case. The choice of the model 
PDE has the purpose to make numerical simulations less time-consuming. On the 
other hand, methods developed in [29, 30, 31] allow naturally direct computation of 
pattern-forming modes without performing any numerical simulation with nonlinear 
PDEs. 
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The knowledge of the structure of resonance clusters might shed some new light 
into the origin of such well-known physical phenomena as Benjamin- Feir instability, 
[3]. As it was shown recently in [18, 19, 51], its explanation as the modulational 
instability, though well established in water waves theory, has to be seriously recon- 
sidered for it over-predicts the growth rate of the waves. Another inconsistency is 
due to the fact that it can easily be stabilized by arbitrary small dissipation. The 
other way to treat the problem would be to explain the modulational instability in 
terms of non-collinear resonances. 

4.1. Quasi-resonances and approximate interactions. The reason why the 
predictions of CWT theory are mostly not corroborated, both in numerical sim- 
ulations and laboratory experiments, is that resonance broadening Q (also called 
resonance width or frequency discrepancy or frequency mismatch etc.), defined as 

Iwi ± W2 ± • • • ± w^l > 0, ki ± ka ± • • • ±k, 0, s < oo, (24) 

is not big enough for CWT theory to be applicable. CWT theory is supposed to 
work when the resonance broadening Q is greater than the spacing S^^ between the 
adjacent wave modes 

n > {dLj/dk)2TT/L, (25) 

where L is the box size. The trouble comes from the notorious small divisor problem 
known from KAM-theory; in CWT theory it appears first in [67], Eq.(2.5.2). The 
CWT theory assumes weak nonlinearity, randomness of phases, infinite-box limit, 
and existence of the inertial interval {ko,ki), where energy input and dissipation 
are balanced. Under these assumptions, the wave system is energy conserving, 
and wave kinetic equations describing the wave spectrum and energy power spectra 
k~°', a > 0, have stationary solutions ([44, 64, 66], etc.) For <C fco, so-called finite- 
size effects take place which arc due to boundary conditions and should be regarded 
separately. For A; 3> fci, dissipation suppresses nonlinear dynamics (see Fig. 2, left). 
In the last decade, this standard view proved to be incomplete while finite-size 




Figure 2. 



effects are well observable within inertial interval; this led first to introducing such 
special types of wave turbulence as frozen turbulence of capillary waves [48] and 
mesoscopic turbulence of gravity waves [65] . Quite recently, a novel two-layer model 
of laminated turbulence has been developed [27]. It allows to explain finite-size 
effects in the same frame: each small divisor leaves a gap in the power spectrum, 
shown as an empty circle in Fig. 2 on the right. The lower boundary for the radius 
TZ of each circle can be computed using Thue-Siegel-Roth theorem [28] for the case 
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of an irrational dispersion function, which covers all water waves (4)-(5). The gaps 
corresponding to the exact and quasi-resonances are shown by yellow squares; their 
dynamics is covered by the dynamical systems of corresponding resonance clusters 
([24], Fig.l). Some gaps, shown by blue diamonds, correspond to the discrete modes 
that do not take part in resonances; they just keep their initial energies ([24], Fig. 3). 

Definition 4.1. Solutions of (24) are called: exact resonances, if f2 = 0; quasi- 
resonances, if 

< 7^ < n<{duj/dk)2Tr/L; 
and approximate interactions, if 

n > {duj/dk)27r/L. 

Exact resonances and quasi-resonances form discrete layer of turbulence, while ap- 
proximate interactions form continuous layer. Correspondingly, the notion of dis- 
crete wave turbulence (DWT) is nowadays used, as opposed to CWT. 

Remark 5. Previous to our theory [27], the standard counterpart of CWT were 
low-dimensioned systems characterized by a small number of modes included. On 
the other hand, DWT is characterized by the clustering itself, and not by the number 
of modes in particular clusters which can be fairly big. 

The minimal resonance broadening, necessary for quasi-resonances to appear, 
has been investigated numerically in [56] both for capillary and gravity water waves. 
Numerical simulations [54, 55, 56] showed that exact and quasi-resonances among 
the gravity waves are observable not on the time scale 0{t/e^) of kinetic theory but 
on the linear time scale. This means that dynamics of discrete layer but not the 
kinetic equation should be used as the base for a short-term forecast. The classical 
notion of wave interactions would be quite inappropriate - even misleading ~ for 
the description of all these important effects. 

5. Summary. In 1981, in the paper with the title "Wave interactions - evolution of 
the idea," [47], O. Phillips wrote prophetical words: "New physics, new mathemat- 
ics and new intuition is required" in order to gain some understanding of finite-size 
effects in wave turbulent systems. Few years later, two - nowadays classical - books 
were published: "Wave interactions and fluid flows," [12], and "Kolmogorov spectra 
of turbulence," [66]. Appearing in 1985, the first one gives comprehensive account of 
the theory and experiments on exact resonances, quasi-resonances and approximate 
interactions - all together named wave interactions - of primary clusters in three- 
and four-wave systems. The second book came to world in 1992, with a beautiful 
statistical theory of wave turbulence which is due to the exclusion of exact reso- 
nances and quasi-resonances from the consideration. The first attempt to fill in the 
gap between these two treatises has been made in 1994, in a small paper "Weakly 
nonlinear theory of finite-size effects in resonators", [24]. Today, some 15 years 
later, we know much more about these effects and can formulate the three most 
important mathematical problems whose solutions would contribute enormously to 
our understanding of nonlinear resonance dynamics of water waves. 

1st problem. Given a dispersion function and a small parameter, by some stan- 
dard technique initial nonlinear PDE with fixed boundary conditions can be trans- 
formed into a set of dynamical systems corresponding to resonance clusters. Each 
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of them is a Poincare normal form. Classification of all normal forms appearing this 
way according to their integrability properties is needed. 

2nd problem. Wave turbulence with continuous spectra is described quite com- 
fortably by the energy power spectra fc~". stationary solutions of kinetic equations: 
the only variable here is the wave length k, while a is a parameter, defined by the 
wave type. Only the fact that the wave system has a conservation law, say its energy 
E, is used, not the value of E. The finite answer is therefore independent on the 
initial conditions. The situation is different for discrete wave turbulence. Energy 
transfer is presently described by exact formulae (20) which are definitely not easy 
to use while initial values of the conservation laws (17) are included explicitly and 
they define the period of energy oscillation of each resonant mode. Some simpli- 
fied description is needed, perhaps also stationary, with ip regarded as a (different) 
constant at each moment of time. 

3rd problem. To describe the transition from discrete to continuous layer of 
turbulence in terms of energy fiow is a very challenging problem, which probably can 
not be solved before a simplified description of the discrete layer will be obtained. 
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